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Abstract 

A purely atomic immigration superprocess with dependent spatial mo- 
tion in the space of tempered measures is constructed as the unique strong 
solution of a stochastic integral equation driven by Poisson processes based 
on the excursion law of a Feller branching diffusion, which generalizes the 
work of Dawson and Li [3]. As an application of the stochastic equation, 
it is proved that the superprocess possesses a local time which is Holder 
continuous of order a for every a < 1/2. We establish two scaling limit 
theorems for the immigration superprocess, from which we derive scaling 
limits for the corresponding local time. 
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1 Introduction 



Let p > and let = (1 + x^Y^I'^ for x G M. We denote by Cp(]R) the set of 

continuous functions on M satisfying |0| < const ■ 0p and denote by Cp(M) the subset 
of Cp(M) consisting of twice continuously differentiable functions with |0'| + |0"| < 
const ■ 0p. Let Mp(]R) denote the space of tempered Borel measures /i on M such that 

(0, /i) := / 0(x)/i((ix) < oo 

for every G Cp(M). Let Mp(M) be the subset of Mp(]R) consisting of purely atomic 
measures. In the case p = 0, we simply write C(M) and M(]R) instead of Co(M) and Mo(M), 
respectively. Let "|| ■ ||" denote the supremum norm. Suppose that /i is a continuously 
differentiable function on M such that both h and h! are square-integrable. Then the 
function 

p{x) = [ hiy - x)hiy)dy, x eR (LI) 

is twice continuously differentiable with bounded derivatives p' and p". We fix a constant 
cr > and a measure m G Mp{R). Let q^UiO) be a Borel function of (i^, a) G Mp(M) x 
M satisfying certain regularity conditions to be specified. A martingale problem for a 
continuous process {Yt : t > 0} in Mp{R) can be formulated in the following way: For 
each G C2(M), 

Mi(0) = (0, Yt) - (0, Fo) - ^P(O) ^ (0", Y,)ds - ^ (g(y;, ■)0, m)ds (1.2) 
is a continuous martingale with quadratic variation process 

(M(0))i= ! {a<p\Y,)ds+ [ ds [ p{x-y)<P\x)<P\y)Ys{dx)Y,{dy). (1.3) 
Jo Jo JR2 

A solution {Yt : t > 0} of the martingale problem is called an immigration superprocess 
with dependent spatial motion (ISDSM); see Ref. [3]. The phrase superprocess with depen- 
dent spatial motion (SDSM) naturally refers to the special case m = 0. Compared with 
the classical super Brownian motion over the real line, the last term in (11. 2p represents an 
immigration factor with interactive immigration rate while the last term in (11.31) comes 
from the dependent spatial motion. In particular, if q^u, a) = q{a) only depends on a G M, 
the immigration becomes non-interactive and the uniqueness of solution of the martingale 
problem can be proved by a duality argument or a conditional log-Laplace functional; see 
Refs. [U [HI [TOl [16]. The uniqueness of solution of the general martingale problem still 
remains a challenging open problem. 

An ISDSM was constructed in Ref. [3] in the special case p = 0. Instead of the mar- 
tingale problem, the authors considered a stochastic integral equation driven by Poisson 
random measures based on the excursion law of a Feller branching diffusion. They showed 
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that there is a unique strong solution of the equation which also solves (11.21) and (11.31) : 
see also Refs. [71 [13] . In this paper we extend the result of Ref . [3] to the general state 
space Mp(M). For the study of the ISDSM, the stochastic equation has several advantages 
over the martingale problem formulation. For instance, the stochastic equation provides 
much more information on the structures of the ISDSM than the martingale problem. As 
pointed out in Ref. [3], from the equation we know immediately that the ISDSM lives 
in the space of purely atomic measures in contrast to the classical super Brownian mo- 
tion; see also Ref. [15]. With that observation, it is natural to ask whether or not the 
corresponding occupation time process 

Zt:= [ YJs, t>0 (1.4) 
Jo 

is absolutely continuous with respect to the Lebesgue measure. We give an answer to 
this question as an application of the stochastic equation. We shall see that {Zf : t > 
0} is really absolutely continuous and its density field {z{b,t) : 6 G M, t > 0} can be 
represented in terms of stochastic integrals of the excursions with respect to the Poisson 
random measures and certain Brownian local times. By this representation we prove that 
{z{b,t) : 6 G M, t > 0} is Holder continuous of order a for every a < 1/2. It seems 
difficult to establish those results from (II. 2p and (II. 3p . The simple derivations of the 
results given here show the efficiency of the stochastic equation in the study of properties 
of the immigration superprocess. As another application of the stochastic equation, we 
prove two scaling limit theorems for {Yt : t > 0}. From those theorems we derive the 
limit theorems for the density field {z{b,t) : 6 G M, t > 0}. 

2 Preliminary results 

Let m G Mp(M) and let (f2, ^,^t,P) be a filtered standard probability space satisfying 
the usual hypotheses. A function r]{-, -, ■) on [0, oo) x R x is said to be simple if it is of 
the form 

oo 

ri{s,x,u) = rio{x,uj)liro,r^]{s) + '^Vi{x,i^)l{n,ri+r]{s), (2.1) 

i=l 

where = tq < ri < r2 < . . . and f]i{-,-) is i^(M) x ^rr^^sasurable. Let <^ be the 
completion with respect to dsm{da)P{duj) of the a-algebra on [0, oo) x M x f2 generated 
by all simple functions. We say a function on [0, oo) x M x f2 is predictable if it is J^- 
measurable. Let L'^{j^) denote the set of all non-negative predictable functions ?7(-, ■, ■) 
on [0, oo) X M X f2 such that 

[ E[{r]{s,-)(pp,m)^]ds < oo, t > 0. (2.2) 
Jo 



3 



Suppose that a > and ?7(-, ■, ■) G L^(^). Given /i G Mp(M.), we consider the following 
martingale problem of a continuous process {Yf : t > 0} in Mp(M): For each G 



Mi(0) = (0, r,) - (0, fi) - -p(0) (0", Y,)ds - (r/(s, O^, m)ds (2.3) 
is a continuous martingale relative to {^t)t>o with quadratic variation process 

{M{<f)))t= [ {a<p\Ys)ds+ [ ds [ p{x-y)<P\x)<P\y)Ys{dx)YMy). (2.4) 

JQ Jo Vr2 

A solution of (12. 3p and (12.41) can be regarded as an generalized immigration superprocess 
with immigration rate given by the two parameter process {t]{s, a) : s > 0, a G M}. 

Proposition 2.1 Suppose that {Yt : t > 0} is a solution of the martingale problem 
given by ( 12.3)) and Ii2.4\) . Let c > be a constant such that \(f>p\ + < c(f>p and let 
C{t) =c2||p||(16+ ||p||t). Then we have 



E[{^p,Y,)]<Gt{<Pp) + -c\\p\ 



Gs{(pp) exp |^c||p||(t - s)^ds 



and 



where 



and 



E 



snp {(f)p,Ysy < Ht{(f)p) + Cit) Hs{(j)p)exp{C{t){t-s)}ds, 

LO<s<i -I 



Gt{(f)p) = {(l)p, 1^) + / mv{s,-)(j)p,m)]ds 



(2.5) 



(2.6) 



Ht{(j)p) = A{<f)p, pf + At [ E[(r/(s,-)0p,m)2]rfs + 16a / E[(0p, 

Jo Jo 

Proof. The first inequality follows by taking the expectations in (12. Sp and applying 
Gronwall's inequality. By (12.31) and the martingale inequality we have 

(?7(s, ■)(f)p,m)ds 



E 



sup {(Pp,Ysf < 4(0p,/i)2 + 4E 



0<s<t 







IpfE 



+ 16 / E 







a{(l>l,Ys) + \\pM,Ys)' 



ds 



< 4(0p,/i)2 + 4t / E[{r]{s,-)(j)p,mf]ds 
Jo 



+ c'\\pft I E\ 



Ysf]ds 



+ 16 ^ (^aE[{<f)p, Ys)] + c'\\p\\E[{<Pp, n)^])^^. 



Therefore, we can use Gronwall's inequality again to get (12.61) . □ 

Clearly, the family of martingales {M^ ((/>)} in (12.31) defines a martingale measure 
M{ds, dx) in the sense of Ref. [H]. The next result follows by standard arguments. 

Proposition 2.2 Suppose that {Yt : t > 0} is a solution of the martingale problem given 
by dl^) and Then for any t > and G C^{R) we have a.s. 

{<f),Yt) = {Pt<P,fi)+ [ {r]{s,-)Pt-s<P,m)ds+ [ [ Pt-s<P{x)M{ds,dx), (2.7) 
^0 Jo Jr 

where {Pt)t>o is the semigroup of the Brownian motion with quadratic variation p{0)dt. 
Consequently, 

E[(0, Yt)] = {Pt<P, fi)+ [ E[(r/(s, ■)Pt-s<P, m)]ds. (2.8) 

Now we can give the construction of a solution of the martingale problem (12. 3p and 
(12. 4p with initial state /i G Mp(M). Let {B{t) : t > 0} be a standard Brownian motion. 
For any initial condition ,^(0) = x > the stochastic differential equation 

d^{t) = ^Jai{t)dB{t), t > 0, (2.9) 

has a unique solution, which is known as a Feller branching diffusion with constant branch- 
ing rate a. In the sequel, we simply call {^(t) : t > 0} a a-hranching diffusion. Let 
W = C([0,oo),M+) and let to{w) = inf{s > : w{s) = 0} for w e W. We denote by 
Wq be the set of paths w E W such that w(t) = w{0) = for every t > to{w). Let 
(^(PVq)) ^i(W^o)) be the natural a-algebras on Wq generated by the coordinate process 
and let be the excursion law of the a-branching diffusion defined in Ref. [3]. Suppose 
that on a complete standard probability space {Q, ^ , P) the following are defined: 

(2. a) a white noise Wids^ dy) on [0, oo) x M based on the Lebesgue measure; 

(2.b) a sequence of independent cr-branching diffusions {^i(t) '■ t > 0} with deterministic 
initial values ^i(O) = ^j, i = 1, 2, ■ ■ ■; 

(2.c) a Poisson random measure Ni{ds,da,du,dw) on [0, oo) x M x [0, oo) x Wq with 
intensity dsm{da)duQ,i^{dw) . 

We assume that {W{ds,dy)}, {^iit)} and {Ni{ds,da,du,dw)} are independent of each 
other. For t > let be the a-algebra generated by all P-nuU sets and the families of 
random variables 

{iy([0,s] X B),ii{s) ■.0<s<t;B e ^(M),z = 1,2,---} (2.10) 

and 

{Ni{J X A) -.J e ^([0, s] X M X [0, oo)); A e ^t-siWo); 0<s<t}. (2.11) 
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It is known that for any (r, a) G [0, oo) x M there is a unique solution {xr,a{t) : t > r} of 
the equation 

x{t) = a+! [ h{y ~ x{s))W{ds,dy), t > r; (2.12) 

Jr Jr 

see Refs. [U [15]. Indeed, each {xr^ait) : t > r} is a Brownian motion with quadratic 
variation p{0)dt. Let {cj : i = 1, 2, ■ ■ ■} C M and assume ^i4>p{(^i) < oo. For 

?7(-, ■, ■) G L^(<?^) we define the purely atomic measure-valued process {Yt : t > 0} by 

yt = Y,Ut)Ka^{t)+ / / / w{t-s)5:,^^^(t^Ni{ds,da,du,dw). (2.13) 

Jo Jr Jo Jwq 



(We make the convention that /^^ = /(a^]-) 

Theorem 2.1 Tie process {Yt : t > 0} deGned by 112. 13\) has a continuous modification 
in Mp (M) which solves the martingale problem given by A2.3\) and A2.4\) . 

Proof. For each integer k > 1 let r]k(t, a) = rjit, a) A {k(j)p{a)). Then t ^ E[(?7fc(t, ■), m)^] 
is a locally bounded function and 



/' 

^0 



E[(|r/(s,-)-%(s,-)|0p,m)V5^O 



as ^ oo for every t > 0. By Theorem 5.1 of Ref. [3], the process {Y^^^ : t > 0} defined 
by 

= 6W^xo.„jt) + / / / / w{t - s)6^^ ^^t)Ni{ds,da,du,dw) 

I 1^, ' JO JR JO J Wo 



|ai|<fc 



has a continuous modification in Mq (M) and hence in (M). Note also that {y/^^ : ^ > 0} 
is increasing mk>l. For n > k > 1 let = ?7„ - r/fc, /ife,n = Efc<|a,|<n^*(0)'^«. ^i^*^ 



= z2 ^i(^)^^o,a.(t) + / / / / w{t- s)6^^^^t)Ni{ds,da,du,dw). 

As in the proofs of Theorems 4.1 and 4.5 in Ref. [3] it follows that {xf'^'"'* : t > 0} has a 
continuous modification in Mp(M) and for each (j) G Cp(M), 

")(0) = - - [\g,^^(s,-)<P,m)ds-lp{0) [\<P",Xi''-^)ds 

Jo ^ Jo 

is a continuous martingale with quadratic variation process 

V0',^f'"V^+ / p(x-y)</.'(x)0'(y)Xf'")(dx)Xf'")(rf|/); 

Jo JR2 
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see also Lemma 3.1 of Ref. [7]. Let C{t) be defined as in Proposition 12.11 We liave 



'"^)] < + ^cIIpII I^Gi''-\^,) exp {ic||p||(t - s)]ds 



and 



E 



sup(0p,Xf'"))2 <Hl''''\<P) + C{t) I Hf^-\4>)eMCm-s)]ds 



0<s<t 



where 



G 



{k,n)i 



and 



(k,n)i 







Clearly, we have G 



{k,n). 



and Hr">{(l)p) ^ as A;, n oo. It follows that 



16(7 / E[< 

r(fc,n) , 



E 



sup 



0<s<t 







as /c, n — >■ oo. Then there is a continuous process {Xt : t > 0} in Mp(M) such that 



E 



sup (0p,X,-FW)^ 

0<s<f 



as n — s> OO for every t > 0. Clearly, {Xj : t > 0} is independent of the particular choice 
of the approximating sequence {77^} C L\{!^) and the martingale characterizations (12.31) 
and ( 12.4p hold with {Yt : t > 0} replaced by {Xt : t > 0}. Then, to finish the proof 



we only need to show a.s. Xt = Yt for every t > 0. Recall that J^^^'"^^ means J, 



(0,r,(s,a)] 



m 



fl2.13p . Then we have a.s. Xt < Yt from the limit procedure for the construction of Xt. 
Let {xj.''^ : t > 0} be constructed by the same procedure for ?7(-, -, ■) + 1/k G L'^{j^). 
Obviously, we have a.s. Yt < xf'\ However, by Proposition 12.21 we have 



E[{<p,Xt)] = {Pt<P,fi)+ [ EMsr)Pt-s 

Jo 



and 



EK0,xf'))] = (P,0,^)+ rE[((r/(s, 

^0 



0, m)]ds 



l/k)Pt^s4>,m)]ds. 



It follows that 

E[(0,X,)] < E[(0,F,)] < lim E[(0,xf )] = E[(0,X,)]. 

fe— »oo 

That proves Xt = Yf a.s. for every t > 0. □ 

A solution of fl2.3p and (12. 4p with general initial state fi G Mp{M.) can be constructed 
in the following way. Let {W{ds,dy)} and {Ni{ds, da, du, dw)} be given as in (2. a) and 
(2.c) and suppose we are also given 

(2.d) a Poisson random measure NQ{da,dw) on M x Wq with intensity fi{da)Q,i^{dw). 

We assume that {W{ds,dy)}, {No{da,dw)} and {Ni{ds, da, du, dw)} are independent of 
each other. For t > let be the a-algebra generated by all P-nuU sets and the families 
of random variables 

{W{[0, s] X B), No{F X A) : < s < t, B e ^{R), A e ^t{Wo)} (2.14) 

and 

{iVi(J X A) : J E ^([0, s] X R X [0, oo)), A E ^t_,(lVo), 0<s<t}. (2.15) 
Let {Xf : t > 0} be defined hj Xq = fj. and 

Xt = [ [ w{t)5,,^^it)No{da,dw), t>0. (2.16) 
Jr J Wo 

For r7(-, -, ■) E L\{^) we define 

Yt = Xt+ [ [ r^ '^ f ^{t- s)5.,^^^^t)N^{ds, da, du, dw), t > 0. (2.17) 
Jo Jr Jo J Wo 

By similar arguments as in the proof of Theorem 12.11 we obtain 

Theorem 2.2 The process {Yf : t > 0} defined by Ii2.17\) has a continuous modification 
in Mp(]R) which solves the martingale problem given by ( 12.3)) and \2.4\i . 



3 Stochastic equations for superprocesses 

In this section, we give the construction of the ISDSM by solving stochastic equations 
driven by Poisson random measures on the space of excursions. For any /i and v E Mp(]R) 
set 



1/^ 



sup 



f{x)(t)p{x)^i{dx) - \ f{x)(^p{x)v{dx 



(3.1) 
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where -Bi(ffi) denotes the set of Borel functions / on M such that < 1 for all x G M. 

Suppose that q{-,-) is a Borel function on Mp(R) x M such that there is a constant K such 
that 

{qifi, m) < Kil + {<j)p, /i)), fiE M(R), (3.2) 

and for each R > there is a constant Lf> > such that 

■) - C[{i^, -Mp^m) < LrWij, - iy\\p (3.3) 

for fi and z/ G Mp(M) satisfying {4>p, fi) < R and {4>p,i^) < R. Given yU G Mp(M), we 
consider the following martingale problem of a continuous process {Yt : t > 0} in Mp(M): 
For each G ^^(M), 

M,(0) = (0, r^) - (0, /i) - lp(0) j\<P\ Y,)ds - j\q{Y,, 00, m)ds (3.4) 
is a continuous martingale with quadratic variation process 

(M(0))i= ! {a<p\Z)ds+ [ ds [ p{x-y)<P\x)<P\y)Ys{dx)Ys{dy). (3.5) 



By similar arguments as in the proof of Proposition 12.11 we have the following 

Proposition 3.1 Suppose that {Yt : t > 0} is a solution of the martingale problem 
given by / I3.4)) and ^3.5\) . Let c > be a constant such that |0p| + |0p| < c0p and let 
Ci = K + c\\p\\/2 andC2{t) = c^HpH (16 + ||p||t). Then we have 

E[(0p, Yt)] < Gti<pp) + Ci [ a(0p) exp {Ci{t - s)}ds (3.6) 

Jo 



and 

ft 



E 



sup (0p, Ysf] < Ht{(j)p) + C2it) f H,{(j)p) exp{C2{t){t - s)}ds, (3.7) 

)<s<t -I Jo 



'-0<s<t 

where Gt{(t>p) = (0p, fi) + Kt and 



Ht{<l)p) = 4{<Pp,fiy + AKH' + 8{Kt + 2a) [ E[{<Pp,Ys)]ds. 

Jo 

Let {W{ds,dy)}, {^j(t)} and {Ni{ds,da,du,dw)} be given as in (2. a), (2.b) and (2.c). 
Let {^t)t>o be the filtration generated by fl2.10p and (12. lip . By a solution of the stochastic 
equation 

Yt = y2^i(^)^^o.a,(t)+ / / / w{t- s)5^^jt)Ni{ds,da,du,dw), (3.8) 

j^]^ Jo JrJo JWn 



we mean an (^t)-adapted continuous process {Yt : t > 0} in (M) that satisfies (13. Sp 
for all t > 0. 
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Theorem 3.1 There is an unique solution {Yt : t > 0} of Ii3.8\) . which is a strong Markov 
process in (M) and solves the martingale problem given by A3.4\) and Ii3.5\) . 

Proof. Suppose that {Yt : t > 0} and {Y"/ : t > 0} are two continuous solutions of 
(K8^ . Fix > 1 and let r = inf{t > : F*) > i? or (0p, F/) > R}. Observe that 
q{Yt,x)l{t<T} and x)1{(<t-} are predictable and define 

ptAT p pq{Ys,a)Aq{Y^,a) p 

X* = / w(t - s)6a;,^,(t)Ni{ds,da,du,dw) 

Jo Jr Jo J Wo 

and 

ptAT P pq{Ys,a)Vq(Y^,a) p 

Y* = / w{t - s)5^^^^{t)Ni{ds,da,du,dw). 

Jo Jr Jo Jwo 

By Proposition 12.21 we have 

E[(0,X;)]= / E[{{qiY„-)AqiY:,-))Pt-s<P,m)^s<r}]ds (3.9) 

and 

E[(0,F;)] = f E[{{q{Y„-)y q{Y:,-))Pt-s<P,m)^,<^y]ds. (3.10) 
Jo 

Let c > be a constant such that < 2c(f)p. Then we have 

f^PtM^) = lpt(i>;ix) < cPtMx). 

By a comparison theorem, we get Pt(f)p < e'^^(f)p for every t > 0. Observe also that 
- F^rllp < {(l)p,Y;) - Then and (^JO^ imply that 

nWYL-ytArllp] = f E[{\q{Y:,-) - q{Ysr)\Pt-s<Pp,m)l{,<^y]ds 
Jo 

<e^* / E[{\q{Y:,-) - q{Ys,-)\<Pp,m)^s<r}]ds 
Jo 

<LRe''! E[\\Y'^^^-YsAr\\p]ds. (3.11) 
Jo 

By Gronwall's inequality we conclude E[||Y'/^^ — I^atIIp] = 0. Since i? > can be arbitrary, 
that gives the uniqueness of solution. To show the existence of a solution, we first assume 
(13. 3p holds with replaced by a universal constant L independent of R. We define the 
sequence continuous processes {Y^^^ : i > 0} by setting Y^;^^ = 6(^)'^a:o,ai(t) 

yW^yW^/ / / / nj(t-s)S^^^^^t)Ni{ds,da,du,dw) (3.12) 

Jo Jr Jo JWq 
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for n > 1. By a reasoning as in fl3.1ip we see that 

E[||yW _ < Le^t /"* E[||Fi"-^) - Y}--'X]ds. 

Jo 

By a standard argument, one sees there is a predictable process {Yt ■ t > 0} such that 
hm^^oo E[||yf — y}"'^ lip] = uniformly on each bounded interval of t > 0. By Theorem 12 .11 
we see that 

Yt = Y}''^+[ [ r^^'^ f yj{t-s)6^^^^^t)Ni{ds,da,du,dw) (3.13) 

Jo JrJo J Wo 

defines a continuous process {Yt : t > 0} on Mp(M). Based on fl3.12p and (13.131) we may 
follow the calculations in (13. lip to obtain 

E[r/"^ - YtU < Le^' fnWYj--'^ - YX]ds. 

Jo 

Letting n — ^ oo we see that Yf = Yt a.s. for every t > 0. That proves the existence 
of a solution of (13.80 for a universal constant L. The extension of the existence to the 
general condition (13.31) is a standard localization argument. The strong Markov property 
of {Yt : t > 0} follows from the uniqueness of the solution. □ 

To consider a general initial value /i G Mp(M) let {W{ds,dy)}, {NQ{da,dw)} and 
{Ni{ds,da, du,dw)} be given as in (2. a), (2.c) and (2.d). Let {.^t)t>o be the filtration 
generated by fl2lD and (12775]) . Let {Xt : t > 0} be given by fl2T6D . By a solution of the 
stochastic equation 

ft r rq{Ys,a) r 

Yt = Xt+ / wit- s)5.,^^^(t)Ni{ds,da,du,dw), t > 0, (3.14) 

Jo Jr Jo J Wo 

we mean an (^t)-adapted continuous process {Yt : t > 0} in Mp(M) that satisfies (13.140 . 
By arguments similar to those in the proof of Theorem 13.11 we have 

Theorem 3.2 There is a unique solution {Yt : t > 0} of 113. 14\) . which also solves the 
martingale problem given by Ii3.4\) and / I3.5)) . 

4 Existence and continuity of local times 

In this section, we prove that the occupation time of an immigration superprocess {Yt : 
t > 0} has a Holder continuous density field {z{b,t) : b E M.,t > 0}. As in the classical 
case the two parameter process {z{b,t) : 6 G M, t > 0} can be interpreted as the local 
time of {Yt : t > 0}. We shall need the following two lemmas. 
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Lemma 4.1 Let {Bj{-) : j = l,---,n} be a family of Brownian motions, {^j(-) : j = 
1, ■ ■ ■ , n} be a family of independent a-branching diffusions and {<yj{-) : j = 1, - ■ ■ ,n} be 
a family of bounded processes. Suppose that the two families {-Bj(-) : j = I, ■ ■ ■ ,n} and 
{O(') • i ~ 1; ■ ■ ■ ; ''^j ^-^6 independent of each other. Let denote the local time of 

Bj{-). Then to each integer k > 1 there corresponds a constant Ck>0 such that 



j=l -^rj 



" rt 



(4.1) 



for any intervals [rj, tj] C [r,t], j = 1, ■ ■ ■ ,n. 



Proof. In this and the following proofs, Ck will denote positive constants that may change 
values from line to line. By integration by parts formula, we have 



^j{s)aj{s)lj{b,ds) = ^j{tj) I aj{s)lj{b, ds) 

d^jis) I aj{s)lj{b,ds). 
By Holder's inequality and the boundedness of {ttj(-) : j = 1, ■ ■ ■ , n}, 

3=1 ■'■'^ 



(4.2) 



2k-l 



2k 



i=i i=i 
<C,(t-r)'=E|[^e,(t,)J }, 



(4.3) 



where we also used the estimate E[/j(6, [rj,tj])^'^] < Ckitj — and the independence of 
and /j(-, ■) for the last inequahty. Similarly, by Doob's martingale inequality we get 



E 



IfE t di,{s) r a,{s)l,{b,ds)Y'} 



j=l -^'J 

n „f 

<c.e{[x:/' 

j=l -Jrj 



aj{s)lj{b, ds) j ^j{s)ds 
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r j 

<C,e{[J] / l,{h,[r,,s\)%{s)d^ 

i=i "^'■^ i=i "^""^ 

1=1 J-<-i 



Then the desired estimate follows from (14.21) . (14.31) and (14. 4p . □ 

Lemma 4.2 Under the assumptions of Lemma \4.1\ to each integer k > 1 there corre- 
sponds a constant > such that 

^j{s)aj{s)lj{bi,ds) - ^j{s)aj{s)lj{b2,ds) | 

j = l -^^3 -^^3 

<Cu{h,-h,f[Y.[(^Y.^,{t,)) ]+e[(5^/ (4.5) 



for any 6i, 62 G arid arij intervals [r^, tj] C [r, t], j = 1, ■ ■ ■ , ra. 

Proof. By a formula at page 211 of Ref. [12] it is easy to show that 

E{[/,(6i, [r,,t,]) - /,(62, [r„t,])]''} < C,{h - hf . 

We can apply (14. 2 p to the two integrals on the left hand side of (14.51) . Then the result is 
obtained by similar estimates as in the proof of Lemma 14. 1[ □ 

Theorem 4.1 Let ri{-,-,-) G L'^{^) and suppose there is a deterministic increasing 
function f/(-) ori [0, 00) such that r](t, a, u) < fj{t) for all {t, a, u). Let {Zt : t > 0} denote 
the occupation time of the process {Yt : t > 0} defined by 112. 17\) . Then Zt is a.s. absolutely 
continuous and the corresponding local time is given by 

z{b,t) = / / No{da,dw) / w{u)lo^aib, du) 

Jr J Wo Jo 

r-t r- p{s,a) r- r-t 

+ / Ni{ds, da, du, dw) / w{v — s)ls^a{b, dv), (4.6) 

Jo Jr Jo J Wo Js 

where {ls,a{b,u) : u > s} is the local time of {xs^a{u) : u > s} at b & M.. Moreover, the 
two parameter process {z{b, t) : 6 G M, t > 0} has a version which is Holder continuous of 
order a for every a < 1/2. 
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Proof. The existence of the local time follows by (12.171) and Fubini's theorem. Let zo{b, t) 
and Zi {b, t) denote respectively the first and the second terms on the right hand side of 
(14. 6p . We shall only give the proof of the continuity result for zi{-,-) since the proof for 
zo{-,-) is similar. 

Step 1. Let us consider the special case with ji and m G M(]R). Observe that zi{h, 0) = 
for every 6 G M. For any t > r > and 6 G M we have 

ft f i'ri(s,a) r rt 



zi{h,t) - zi{b,r) 




Ni{ds,da,du,dw) / w{v — s)ls^a{b, dv) 



+ 




Wo 

ri{s,a) f r-t 

Ni{ds, da, du, dw) / w{v — s)ls^aib, dv). 

Wo Jr 



For e > let W, = {w e Wq : w{e) > 0}. Then A^i([0, t] xMx [0, r/(t)] x W,) is a.s. finite for 
every t > 0. Let S{q, Ni) denote the intersection of the support of Ni{ds, da, du, dw) with 
the set {{s,a,u,w) :s>0,aG]R, 0<n< f]{s),w G W^}. We can enumerate S{q,Ni) 
into a sequence {{sj,aj,Uj,Wj) : j = 1,2,- ■■} so that < Si < S2 < ■ ■ ■. As in the 
proof of Lemma 3.2 in Ref. [5] one sees that, given {{sj, aj, Uj, Wj{e)) : j = 1, 2, ■ ■ ■} each 
{wj{u) : M > e} is a a-branching diffusion independent of the white noise {W {ds , dy)} . 
For any integer A; > 1 we may use Lemma 14.11 and Fatou's lemma to see that 



E 




ri{s,a) 



< lim inf E 



f f 

' Ni{ds,da,du,dw) / w{v — s)ls^a{b,dv) > 
Wo Js ' J J 



+E 



t -I 2k 

^ji^' - Sj)ls„a,{p, dv) 

S ' 

r<Sj <t 
Wj{v — Sj)dv 



< Cfc(t-r)Miminf |e[(^ ^ Wj{t 



r<Sj <t 

<C,(t-r)^'{E 
+E 




nit) 



Wo 



N 2fc' 

w{t — s)Ni{ds, da, du, dw) j 




Wo Jo 



v{t) rt ^ k 

Ni{ds, da, du, dw) / w{v — s)dv 



<C,(t-r)'=|E[(l,Fi)2'=]+E[(l,Zi)^- 

where {Yt\ is the process defined by f l2.17p from {^(t)}, and {Zt} is the corresponding 
occupation time process. By a similar reasoning as the above we have 

''r}{s,a) f nt 2k 



E 




Wo 



Ni{ds,da,du,dw) / w{u — s)ls,a{b,du) 



< lim inf E 



E 

0<Si<r 



Wj{u- Sj)ls^,a,{b,du) I 
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< (7fe(t-r)Miminf |e[(^ ^ Wj{t 

0<s,< 

<c.«-.)'{e[(/7j 



— Sj)du 

m 



+ E 




w(t — s)Ni{ds, da, du, dw) 



Ni{ds, da, du, dw) / w{v — s)dv 

Wo Jr 



<Ck{t-r)'[E[{l,Y,)''] + E[{l,Z,)']y 



Thus we have 



E[|zi(6, t) - z,{b, r)n < C,{t - E[(l, Y.^'] + E[(l, Z,)'] \. 



,2fei 



(4.7) 



For t > and fci,fc2 £ IR we can use Lemma [4.21 and similar arguments as the above to 
show 



E[|zi(fei,t) - zi(62,t)r'] < Ckih - 62)'{e[(1,F<)2^] + E[{l,Zty 



(4. 



Since fi and m G M(M) are finite measures and rjit) is a locally bounded function of t > 0, 
it is not hard to show that E[(l, Ys)'^'' + (1, Z^)''] is locally bounded in t > 0. Then fl47D 
and (14. 8 P imply that Zi{-, ■) has a Holder continuous version of order a for every a < 1/2; 
see e.g. page 273 of Ref. [H]. 

Step 2. Now we consider the general case with /i and m G Mp(M). Let L and T be 
fixed positive constants. For any integer tt, > 1 let a„ > be sufficiently large so that 



I sup |xo,o(s)| > flnj < -• 



(4.9) 



Let 6„ = L + a„ and c„ = I/ + 2a„. Let A„ = {there exists < s < T such that xo,6„(s) = L 
or = Cn} and Bn = {there exists < s < T such that xo,-6„(s) = —L or = — c„}. As 
observed in Ref. [16], any two solutions of (I2.12p started from different locations never 
collide. Then on the event {An U BnY, for any \a\ > Cn and 0<s<M<Twe have 
Xs,a{u) > xo,bn{u) > L OT Xs,a{u) < xq -6„ (m) < -L, and so ls,a{p,u) = whenever \b\ < L. 
It follows that 



z{b,t) 



[-c»i,c„] J Wo 
t 



No{da,dw) / w{u)lo^a{b, du) 



Wo 



Ni{ds,da,du,dw) / w{v — s)ls^a{b,dv) 



for (6, t) G [— -^v, L] X [0, T] on (A„Ui?„)'^. By Step 1, for any a < 1/2 we have a modification 
of z{b,t) on (y4„ U B„y that is Holder continuous of order a in (6, t) G {—L,L) x [0,T). 



15 



In view of (14.91) we have P(A„) = P(-Bn) = 1/n. Then we can modify z{h, t) on the whole 
space so the process becomes Holder continuous of order a in {h,t) G {—L,L) x [0,T). 
Since L > and T > can be arbitrary, we have the desired result. □ 



Theorem 4.2 Let q{-,-) be a Borel function on Mp(]R) x R satisfying the conditions in 
the last section. In stead of / I3.2)) . we assume the stronger condition 

g(/i,a) < J^(l + (0p,/i)2), a e M,/i e Mp(R). (4.10) 

Let {Zt : t > 0} denote the occupation time of the ISDSM defined by 113. 14\) . Then Zt is 
a.s. absolutely continuous and the corresponding local time is given by 

z{b,t) = No{da,dw) / w{u)lo^aib, du) 

Jm. J Wo Jo 



+ / / / / Ni{ds,da,du,dw) / w{v — s)ls^a{b, dv). (4.11) 

Jo Jr Jo J Wo Js 

Moreover, the two parameter process {z{b, t) : 6 G M, t > 0} has a version which is Holder 
continuous of order a for every a < 1/2. 

Proof. Let {y}""^ : t > 0} be defined by the right hand side of (I3.13P with q{Ys, a) 
replaced by ?7n(s, a) := nAq(Ys, a). From Theorem 14. II we know that the occupation time 
of {f/"^ : t > 0} is absolutely continuous with density given by 

Znib,t) = No{da,dw) / w{u)lo,aib, du) 

Jr J Wo Jo 

fqri{s,a) f ft 




Ni{ds, da, du, dw) / w{v — s)ls,aib, dv) 

JRJO JWo Js 

which has a Holder continuous version of order a. The same conclusion is clearly true for 
the occupation time of [Yt : t > 0} on the event {supQ<g<j qiXs, a) < n}. By Chebyshev's 
inequality and the assumption on g(-, ■) is easy to show that 



P| sup qiYs,a) > n| < — e[i + sup 

n<s<t J n I o<s<t 



By Proposition 13. 11 the right hand side tends to zero as n ^ oo. Then we have the desired 
result. □ 



5 Scaling limit theorems 

Scaling limit theorems of SDSM without immigration were investigated in Refs. [5], [9]. A 
direct generalization of the limit theorem of Ref. [5] to the ISDSM was given in Ref. [6]. 
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In this section, we prove two scaling limit theorems for the ISDSM. The limit processes 
obtained here are different from those in the previous work. As consequences, we also 
obtain scaling limits for the corresponding local time. For simplicity we assume 10 = 
and focus on the influence of the immigration. 

Let M = M U {00} be the one-point compactification of the real line and let M(M) be 
the space of Borel measures on M. We fix a metric on M(]R) compatible with the weak 
convergence and regard M (M) as a subspace of M(M) comprising measures supported by 
M. Then a metric on Mp(M) can be defined through the isomorphism 

$p : fi{dx) (-^ (t>p{x)fi{dx) (5.1) 

between Mp(M) and M(M). Let C([0, 00), Mp{M.)) be the space of continuous paths from 
[0, 00) to Mp{M.) endowed with the topology of locally uniform convergence. 

We first assume p > 1 so the Lebesgue measure A belongs to Mp(M). Let q{-,-) be 
a bounded Borel function on Mp(R) x M satisfying the local Lipschitz condition (13.31) . 
Suppose that {W{ds,dy)} and {Ni{ds,da, du,dw)} are given as in (2. a) and (2.d) with 
m replaced by A. Let {xr^ait) t > r} he defined by (I2.12p and let {Yt : t > 0} be the 
solution of 

i-t r- l'qiYs,a) r- 

Yt= / w{t - s)6^^^^^t)Ni{ds,da,du,dw). (5.2) 

Jo JmJo J Wo 

For any integer k > 1 let Y^{dx) = k~'^Yk2i{kdx) . 

Theorem 5.1 Suppose that q{v, a) (l{oo) as \a\ 00 for all fi G Mp(M). Then, ask ^ 
00, {k-^Yl" : t > 0} converges to {q{oo)tX : t > 0} in probability on C([0, 00), Mp(M)). 

Proof. In this and the following proofs, we write "=d" for the equivalence in distribution 
of two processes. We shall also consider some new Poisson random measures which might 
be defined on some extensions of the original probability space. 
Step 1. Let e C2(M). From (I52D we have 

f-t p i-q{Y^2,M) r 

{(P^Yj^) = II I Wk{t- s)(f){x^^f^{t))Ni{Pds,kdb,du,dw), (5.3) 

Jo Jr Jo J Wo 

where Wk{t — s) = k~'^w{k'^{t — s)) and {x'^f^{t) : t > s} is the unique solution of 

x{t)=b+ [ [ hk{y - x{u))k-^/^W{k^du,kdy), t>s 
J s Jr 

with hk{z) = \/kh{kz). It is easy to see that k~^^'^W{k'^du, kdy) is a white noise based on 
the Lebesgue measure and Ni{k'^ds, kdb, du, dw) is Poisson random measure with intensity 
k^ dsdhduQli^{dw) . By the scaling property of the a-branching diffusion it is easy to check 
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that {wk{t) : t >0} under fc^Q^ has the same law as : t > 0} under Q^. From (15.31) 

we get 




JM. JO 



w{t — s)(j){x^ f,{t)) Nk{ds , db, du, dw), 



(5.4) 



Wo 



where Nk{ds,db,du,dw) is Poisson random measure with intensity kdsdhduQ^i^dw) . Let 
ipk = k~^(j). Since QK[w(t)] = 1 for every t > 0, from (15.41) we have 



E[(V^fc,F,'=)] = E I ds db q{Yk2,,kb)w{t-s)(j){xl,{t))Q,{dw) 

Jr Jwo 



ds / E[g(n2„A;6)0«,(t))]rf6 
Jr 

ds db E[q{Yk2s,kb)\x';^^{t) = z](j){z)gp^o)(t-s){b,z)dz, 



where gu(b,z) denotes the density of the heat kerneL By dominated convergence we see 
that 

Mm E[{iljk,Y^^)] = / ds db q{oo)(j){z)gp(o){t-s){b, z)dz = tq{oo){(l), X) . 
Jo Jr Jr 

For every t > we have Qk[w^(^)^] = cr^; see e.g. Lemma 3.1 of Ref. [13]. It follows that 



E 



k-^E 



ds db q{Yk2s, kb)w{t - sy(j){xl,{t)yq^{dw) 

Ir J Wo 



< A;-V||g|| [ {t-s)ds [ E[(j){x';f,{t)y]db 

Jo Jr 

< A;~V||g|| / {t-s)ds / (j){bfdb, 

Jo Jr 



which tends to zero as A; — oo. That proves limfc^oo('^A:7 ^t^) = 9(00)^(0, A) in probability. 
Step 2. From (15. 4p and Theorem 12.11 we see that 



Mf(0) := (0,F,'=)-lp(O) f{<p"X)ds-k [\q{Y,2„k-)<P,X)ds 

^0 Jo 

is a continuous martingale with quadratic variation process 



(5.5) 



{M\<P)), = / {a<P^Y^')ds+ / ds / p{k{x-y))<P'{x)<P'iy)Y^'{dx)Y^\dy). (5.6) 
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Let {rfc} be a bounded family of stopping times. By (15.51) and (15. 6p we have 

+ 3kH [ E[(g(n2(,^+,),fc-)0,A)2]rfs, 
Jo 

where 

Jo Jo 
By Proposition 12.11 it is easy to show that 

linisupE[|(V;fc,Fi+,) - 0. 

k>l 

Another apphcation of Proposition 12.11 shows 

hm snpP{{ipk,Y^'') > a} = 0. 
fc>i 

Then a criterion of Ref. [1] imphes that the sequence {{i/jk,yt^) : t > 0} is tight in 
C([0,cx)),M). Let $p be defined by ([51]). By Theorem 3.7.1 of Ref. [2], {fc-^$pi;^ : t > 0} 
is a tight sequence in C([0, cxd), M(M)). Now the result of the first step implies that 
{k^^^pY^'' : t > 0} converges to {t^pX : t > 0} in probability in C([0, oo), M(R)). 
Since all the processes live in C([0, oo), M(]R)), the convergence also holds in this smaller 
space. In other words, {k~^Y^'' : t > 0} converges to {q{oo)tX : t > 0} in probability in 
C([0,oo),Mp(M)). □ 

The above theorem implies the following scaling limit theorem for the local time of 
the ISDSM. 



Corollary 5.1 Let z{-, ■) denote the local time of {Yt : t > 0} given by the second term on 
the right hand side of A4.11\) . Under the condition of Theorem \5.1\ Zk{t, ■) := k~^z{k -, k'^t) 
converges weakly to g(oo)t^/2 in probability as k oo. 

For a finite reference immigration measure, we can prove a limit theorem which gives an 
interesting random limit process. To describe the limit process we introduce the following 
concept: A two-parameter process {Vrit) : t > r > 0} is called a restricted coalescing 
Brownian flow (RCBM flow) with speed p > provided 

(5. a) for any fixed r > 0, the process {yr{t) : t > r} is a Brownian motion with speed p 
started from yr{r) = 0; 

(5.b) for any fixed s > r > 0, the process {ys{t) — yr{t) '■ t > s} is a Brownian motion 
with speed 2p stopped at zero. 
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Let m G M(M) and let ■) be a bounded Borel function on M(M) x M satisfying the 
local Lipschitz condition (13.31) for p = 0. In addition, we assume 

(5.c) there is a constant e > so that g(i^, a) > e for all u G M(M) and a G M; 
(5.d) g(i/, a) — > q{a) as (1, i^) ^ oo for all a G M. 

Suppose that {W{ds,dy)} and {Ni{ds,da, du,dw)} are given as in (2. a) and (2.d) and 
let {Yt:t> 0} be defined by Let Y^''{dx) = k-^Yk2t{kdx). 

Theorem 5.2 Under the above conditions, as k ^ oo, the sequence {Y^'^ '■ t > 0} 
converges in distribution on C([0, cxd), M(M)) to 

Yj°° := / / / / w{t- s)6^oo^t)Ni{ds,da,du,dw), t>0, (5.7) 

Jo Jr Jo J Wo 

where {x^(t) : t > s} is an RCBM How independent of {Ni{ds,da,du,dw)} with speed 
p(0). 

Proof. Let {xj^(t) : t > s} be defined as in the proof of Theorem 15. 1[ For any (f) G C^(]R) 
we have 

{<P,Y,')= // / Wk{t-s)<P{xl,/,{t))N^{k'ds,da,du,dw). (5.8) 

Jo Jr Jo J Wo 

It is easy to see that Ni{k'^ds,da,du,dw) is a Poisson random measure with intensity 
k'^dsm{da)duQ,K{dw). Let be defined by the right hand side of (15. 2p with q(Ys,a) 
replaced by q{a). Then we have 

{<P,X^)=f f r^'^^ f w,{t-s)<P{xl^/,{t))N^{k^ds,da,du,dw). (5.9) 
Jo Jr Jo Jwo 

The scaling property of the a-branching diffusion implies that 

(0,Xf)=d/ / / / wit- s)(t){xl,,k{t))Ni{ds,da,du,dw). 
Jo Jr Jo Jwo 

Therefore the distribution of : t > 0} on C([0, cxd),]R_|_) is independent oik>l. 

Indeed, by Theorem 4.1 of Ref. [H] it is easy to see that : t > 0} is distributed 

as the cr-branching diffusion with immigration that solves the stochastic equation 

dzt = y/aztdBt + {q, m)dt, zq = 0. 

Given (si, ■ ■ ■ , s^; ^i, ■ ■ ■ , cim), it can be proved as Theorem 2.2 of Ref. [S] that {{x^^ ai/k(')^ 
■ ■ ■ ' -^sm a™/fc('))} converges in distribution to {(x^(-), ■ ■ ■ , (■))}. By dominated conver- 
gence it is easy to show that {(0, Xj^) : t > 0} converges in finite dimensional distributions 
to 

pt p pq{a) p 

(0,r~)= / / / / w{t- s)(f){x^{t))Ni{ds,da,du,dw), t>0. 
Jo Jr Jo Jwk 
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On the other hand, from (15. 8p and (15.91) it is easy to see that 



\{<P,X^)-{<PX)\< / / / Wk{t-s)U\\N,{kHs,da,du,dw) 

Jo Jr Jq{a)Aq{Y^.2^,a) JWq 

and so 

E[|(0,Xf) - {<P,Y,')\] < II0II fds [ E[\q{a)-q{Y,2,,a)\]m{da), (5.10) 

Jo Jr 

By Condition (5.c) we find that (1, 1^2^) = k'^{l, Y^) ^ oo in probabihty for every s > 0. 
Then (15.101) and Condition (5.d) imphes that \{(j),X^) — {(j),Y^'')\ ^ in probabihty. 
Therefore {{(f),Y^'') : t > 0} also converges to {{(f),Yf°) : t > 0} in finite dimensional 
distributions. The tightness of {F/' : t > 0} can be established as in the proof of Theo- 
rem [5]T1 so the sequence converges to {Y^°° : t > 0} in distribution on C([0, oo), M(M)). 
□ 

Corollary 5.2 Let z{-, ■) denote the local time of {Yt : t > 0}. Under the condition of 
Theorem \5.2\ k~^z{k-,kH) converges weakly in distribution to 

z°^{-,t):= f f r'^^ f Ni{ds,da,du,dw) [ w{v - s)l^{-,dv), (5.11) 
Jo Jr Jo Jwq Js 

where {l'^{b,u) : u > s} denotes the local time of {xf{u) : m > s} at 6 G M. Moreover, 
the process {z°°{b,t) : 6 G M, t > 0} has a version which is Holder continuous of order a 
for every a < 1/2. 

Proof. The first part is an immediate consequence of Theorem 15. 2[ The second part 
follows by arguments similar to those given in the proof of Theorem 14. 1[ □ 
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